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Abstract 

In the case of minimal supersymmetric extension of the Stan- 
dard Model (MSSM), when the pseudoscalar Higgs boson mass is less 
than the supersymmetry energy scale, the effective theory at the elec- 
troweak scale is a two-Higgs-doublet model. We diagonalize the mass 
matrix of the general two-Higgs-doublet model, expressing Higgs bo- 
son self-couplings in terms of two mixing angles and four Higgs boson 
masses, and derive in a compact form the complete set of Feynman 
rules, including quartic couplings in the Higgs sector, for the case of 
CP-violating potential. Some processes of double and triple Higgs 
boson production at a high-energy linear collider are calculated in the 
case of mixing angles and scalar boson masses satisfying the MSSM 
constraints. 
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1 Introduction 



A particularly simple extension of the Standard Model containing two scalar 
doublets P] has been very extensively investigated in the framework of min- 
imal supersymmetry. In order to cancel gauge anomalies introduced by the 
fermionic superpartners of gauge bosons and to generate masses of up- and 
down- quarks in a consistent manner two doublets of Higgs fields are neces- 
sary. 

Soft supersymmetry-breaking terms introduce large radiative correc- 
tions to the tree-level Higgs boson masses and couplings H and the effective 
lagrangian of the Higgs sector at the electroweak scale does not satisfy the 
supersymmetry constraints valid at the SUSY scale. In the most general 
case when the supersymmetry scale and the scale of heavy Higgs boson mass 
(usually defined by the mass of the pseudoscalar) are different, the effective 
theory at the electroweak scale is a two-Higgs-doublet model, where the self- 
interaction couplings are defined by the renormalization group evolution of 
the supersymmetric potential couplings from the SUSY scale down to the 
electroweak scale ||, ^ . 

The investigation of direct phenomenological consequences of a two-doublet 
Higgs sector at a future high luminosity colliders, such as LHC and TESLA, 
could provide a possibility to study in detail the structure of effective Higgs 
potential, mass spectrum and couplings of the scalar particles. As usual, 
the variety of channels where scalars could be produced individually or in 
association with vector bosons requires a systematical calculation in order to 
find out what particular channels could have a sufficient counting rates for 
experimental detection at a given collider luminosity. 

We propose a convenient compact form of Feynman rules for a general 
two-Higgs-doublet model that can be used in the following systematical study 
of the Higgs boson production channels and use these rules for the calculation 
of two and three Higgs boson production at a high energy e + e~ collider. 

2 Diagonalisation of the mass matrix in the 
general two-Higgs-doublet model 

General form of the (nonsupersymmetric) SU(2) x U(l) invariant potential 
in the case of two doublets of complex scalar fields <pi, y9 2 can be found in || 

v(vi, ¥*) = Mvin - i f + Writ* - f ) 2 (!) 

+A 3 [Kv?i-f) + (^ 2 -f)] 2 
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+A 4 [OM)(>2>2) - {<Pl<P2)(<P2<Pl)] 

+X 5 [Re((ft<f 2 ) - ^Re(e^)} 2 + A 6 [/m(<^ 2 ) - f/m(e^)f 
where \ are real constants. Components of scalar doublets (pi }2 are 

lfl = {-i w +^^-( Vl J r h l +izi)}, ipi = {-iwf,-^=(v 2 + h 2 + iz 2 )}. (2) 

where if is a complex field and z, h\ y2 are real scalar fields. Vacuum expec- 
tation values t>i, v 2 correspond to the minimum of the potential 

Pi = -^{(W, ^ 2 = -J={0,i; 2 e 4 «} (3) 

where the phase £ can be removed by the rotation of <fff2 not affecting 
the A4 term in (1). Substiution of (2) to (1) gives a bilinear form of the 
mass term with mixed components w, hi t2 , z, which can be diagonalized by 
an orthogonal transformation of the fields in order to define the tree level 
masses of physical bosons. The resulting spectrum of scalars consists of two 
charged H ± , three neutral h, H, A scalar fields, and three Goldstone bosons 
G. This procedure is described in many papers (for instance, f| |7|]). The 
Wi : 2 sector is diagonalized by the rotation of W\, w 2 — > H, G 

wf = -H^sp + G ± c /3 , wf = E ± c fi + G ± s f3 (4) 

defined by the angle 

(5) 

and leading to the massless G field and the field of massive charged Higgs 
boson H ± , m 2 H ± = \±(vf + v%)/2. The z\, 2 sector is diagonalized by the 
rotation Zi,z 2 — > A°,G defined by the angle (3 (4) and giving again one 
massless field G and the field of CP-odd Higgs boson A with the mass 
m A = ^5( v i + Finally, the hi, h 2 sector is diagonalised by the rotation 

hi, h 2 — > h, H defined by the angle a 

2m 12 mn — m 22 

sin2a = , cos2a = (6) 

'(mn - rn 2 2) 2 + 4m 2 2 y (m n - m 22 ) 2 + 4m 2 2 



where 



m n = i[4u?(Ai + A 3 ) +v 2 X 5 ] 
m 22 = ±[4v 2 (\ 2 + \ 3 ) + v 2 \ 5 ] 
m 12 = i(4A 3 + X 5 )viv 2 
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giving two massive fields of CP-even Higgs bosons H, h with the mass values 



m H,h = m n +m 22 ±\l (mn - m 22 ) 2 + 4mf 2 (7) 

In the explicit form the diagonal mass matrix of scalar fields and the physical 
boson interaction vertices are obtained after the following substitution of \ 
to the potential V(ipi,ip 2 ) (1): 

*i= W\\% c ^ml- f s a _ p m%] + f|A 5 (8) 



4 = -j2 m H± 
6 = i? m A° 



where we used the notation v 2 = v \ + u|, s a = sina, c a = cosa. Diago- 
nalization of the mass term takes place if A5 is arbitrary, but the necessary 
condition for the CP- invar iance of potential (1) is A 5 = A 6 0. Unfortu- 
nately, after the substitution of (8) to the potential (1) the intermediate 
expressions for the four scalar boson interaction vertices turn out to be ex- 
tremely cumberous and it is very difficult to reduce them to some compact 
convenient form, where the dependence of the coupling from the parameters 
could be clearly seen. This is a technical problem of the symbolic manipu- 
laion program || that we used. However, symbolic transformations of the 
intermediate expressions are simpler, if we rewrite the potential V((pi, <p 2 ) i n 
the oftenly used representation 

U ((ft, y 2 ) = -nK^f^) - (4((pt(p 2 ) - fi 2 12 {¥i¥2 + (pt<Pi) (9) 

+X 1 ((pf(p 1 ) 2 + A 2 (<^2 + \ 3 (vtfi)(^ 2 ^2) 
+A 4 0^ 2 )(^2>i) + %[(<Pi<P2)(<Pi<P2) + (vJ^O^aVi)] 

It is easy to check that in the case of zero <pfipi phase the potentials (1) and 
(9) are equivalent if the constants \, [i and Aj are related by the formulas 

Ai = A! + A 3 , A 2 = A 2 + A 3 , A 3 = 2A 3 + A 4 , (10) 

\ a 1 5 . 6 t As A 6 
Aa — — aa + r 1 Ak — — 

2 2 2 2 

and 

^12 = ^5^7?, nl = Xivl + X 3 v 2 + X 3 vl, nl = X 2 v 2 + \ 3 i>l + X 3 vj (11) 
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The expressions (11) are sometimes called 'minimization conditions', if one 
starts from the potential U((pi, ^2) (9), where the symbolic structure does not 
show clearly a possible minimum. In the MSSM A 5 =0 and it follows that //f 2 
is fixed and equal to rn^spCp. If this equality is not satisfied (or, equivalently, 
A 5 7^ A 6 in (1)), CP-violation in the Higgs sector can be introduced. The 
diagonal form of U((px, y9 2 ) and the physical scalar boson interaction vertices 
are obtained by the substitution of the following expressions for Aj and \ii in 
the potential (9): 



Ai = 
A 2 = 


3M(5)X + (£ 


fm\ - 

)H - 


^n 2 1 
0/3 ,i 2 1 


A 3 = 


\[2m 2 H± -J&-- 

V L tl SRCR 


S2/3 V H 




A 4 = 




-2m 2 H± ) 


A 5 = 

i4 = 


1 C **12 


■m\) 


2tgPi4 2 ] 


i4 = 






2ctgf3 /if 2 



Our expressions for A4 and A5 are the same as given in || for the case of 
zero A6 and A7. Complete sets of Feynman rules (unitary gauge) for the 
triple and quartic Higgs boson interactions in the general two-Higgs-doublet 
model with a possibility of CP-violation in the Higgs sector (defined by fi^ 
parameter), are shown in Tables 1-2. These sets were obtained by means of 
LanHEP package M. LanHEP package is a specialized symbolic manipula- 
tion system capable to generate Feynman rules for the SU(2), SU(3) gauge 
invariant lagrangians with arbitrary sets of particle multiplets, in the stan- 
dard input lagrangian format of CompHEP package 0. We do not show 
here a rather long set of Feynman rules in the 'tHooft-Veltman gauge, that 
can be also generated after the introduction of ghost and ghost-goldstone 
lagrangian terms to LanHEP program. [] 

We assume that in the Yukawa sector < (p% > couples only to down 
fermions 

v ud - J md — $i(l + 75)^1 + ^2(1 - l^M] (13) 

2\J 2m w s w cp 

(here for the u, d quarks ipi = {u, V u dd+V us s+V U bb}, ip2 = d and analogous 
structures for s,b quarks and leptons, in the case of quarks V a b denotes the 

1 The generation process takes 15 sec. of CPU time (i686). Complete 
lagrangian tables in CompHEP format and Lan HEP package are available at 
http : //theory . npi .msu. su/~semenov/lanhep .html 
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CKM matrix elements), and < v?2 > couples only to up fermions (so-called 
model of type II 



CTTl — — 

—= [^i(l + 75)^2^2 + ^2(1 - 75)^2^1^2] (14) 

2V ImwSwSp 

(here ipi = {u, V u dd + V us s + V U J)}, ip2 = u and analogous structures for c 
and t quarks). Higgs-gauge boson interaction is defined by the straightfor- 
ward extension of the covariant derivative in the case of two scalar doublets. 
It is easy to find the relation between the vacuum expectation values of the 
potential and the IV-boson mass and coupling g = e/sin^vK 

2 2,2 2 /,r\ 

v = Vl + v 2 = —^ L s w (15) 

following from the structure of scalar fields kinetic term D^ipD^ip. 

From the phenomenological point of view the general multiparametric 
two-Higgs-doublet model is too flexible to be systematically used for data 
analysis. Practically no limits on the masses of individual scalars can be 
set if their couplings to gauge bosons and fermions depend on some free 
parameters, and can be very small in a rather large regions of parameter 



space. Recent discussion of the possible limits can be found in |TTJ . However, 
the parameter space can be strongly restricted by the constraints imposed 
by the supersymmetry. 

Let us consider the reduction of the general two-doublet model Feyn- 
man rules shown in Tables 1,2 to the case of minimal supersymmetry model 
(MSSM). The potential ^(^1,^2) (1) contains eight parameters: two VEV's 
Vi, t>2 and six A, (i=l,...6). Eight parameters of the potential U(<pi,<f2) (9) 
A*i5 A*2; A*i2 an d Aj (2=1,. ..5) can be found using (10), (11). From the other 
side, in order to define the Higgs sector we need eight physical parameters: 
the mixing angle (3 and H^-boson mass mw, mixing angle a, the parameter 
fi\2 and four masses of scalars rrih, m#, tua, m ± . Two VEV's can be ex- 
pressed through mw, tg/3 by (5) and (15) and only one degree of freedom 
remains here. In the case of superpotential five additional constraints are 
imposed, relating all Higgs boson self couplings Aj, (i—l,...5) to the gauge 



coupling constants at the energy scale M SUS y 01 : 



t t g 2 + gf 1 g 2 - g| t g 2 t n n ^ 

Ai = A 2 = — - — , A 3 = — - — , A 4 = - — , A 5 = 0. (16) 

As we already noticed, if A5 =0, //i2 is fixed and CP-parity is conserved. The 
remaining two independent parameters may be used to define all Higgs boson 
masses and mixing angles. One can choose, for instance, r l5 r 2 parametriza- 
tion [13] (r 12 = m \ n I m %) or the well-known m^, tg/3 parametrization. In 
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order to reduce the general two-Higgs-doublet model vertices to the case of 
MSSM it is convenient to use the a, f3 parametrization: 

2 2 0/ o Sin (« + /?) 2 2 n C0S (a + P) 

m h = m z cos2/? — — , m H = m z cos2/3 — , 17 

su(a — p) cos(a — p) 

2 2 sin2(a + (3) 2 2 ■ o a 

m A = m z—^7 m > ^12 = m A sin(3cos(3. 

sinz(a — p) 

Substitution of these expressions to the vertex factors in Tables 1,2 after 
trivial trigonometric transformations reduces them to simpler MSSM factors 
IJ. Complete list of Feynman rules at the MSSM scale is shown in Table 3. 

Renormalization group (RG) evolution of the coupling constants Aj from 
the energy scale Msusy to the electroweak scale Mew violates the constraints 

(16) || and the effective low energy potential at the scale M EW is the po- 
tential of a general two-Higgs doublet model with RG evolved couplings A,. 
At a given values of tua, tg/3 (or a, 0), masses of Higgs bosons and the 
mixing angle a (or itia) at the energy scale Msusy can be obtained using 

(17) . Detailed analysis of the following RG evolution and the calculation of 
leading-logarithmic radiative corrections to the mixing angles, masses and 
couplings of Higgs bosons can be found, for instance, in ||. We briefly 
point out that the additional input parameters to be defined in order to fix 
the scheme are the scale of SUSY breaking Msusy, the mass parameter in 
higgsino-gaugino sector /i, and the squark mixing parameters A. 



3 Multiple production of neutral Higgs bosons 

The processes of multiple neutral Higgs boson production in the MSSM were 
considered in JL4|, |15| in the framework of effective potential approach to the 
calculation of radiatively corrected scalar masses and couplings [|16|] of the 
SUSY Higgs sector. The reactions 

e + e~ — > hhZ, e + e~ — ► hhA, e + e~ — > v e v e hh (18) 

were considered and it was shown that the cross sections of double and triple 
Higgs boson production are not small and the experimental measurements 
of triple Higgs boson couplings are realistic. 

We used the results of ]5], [Tj| to calculate the radiatively corrected masses 
of Higgs bosons and mixing angle a in the renormalization group approach to 
the Higgs potential couplings evolution from the SUSY scale M S usy down 
to the electroweak scale (see also jl7|). We set the Msusy =1 TeV and 
have not included the effects of squark mixing by setting the parameters 
A and /x equal to zero. In the case of not too large tg/3 (we used tg/3 = 
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3) and the pseudoscalar mass rriA of order 150-250 GeV, masses of heavy 
CP-even Higgs boson H and charged Higgs boson H ± are also at the scale 
150-250 GeV. The lightest Higgs boson mass is approaching 100 GeV when 
the pseudoscalar mass surpasses 200 GeV. (Changes of the SUSY scale and 
mixing parameters can in principle shift rrih by about 50 GeV, see the details 
in |5], |15|). These radiatively corrected parameters were used in our set of 
Feynman rules. The following calculation of the complete tree level amplitude 
for the multiple Higgs boson production processes (18) was performed by 
means of CompHEP package ||, when the exact symbolic result for the 
matrix element squared is converted to FORTRAN code and integrated by 
multichannel Monte-Carlo method. The s-channel resonant peaks of the 
amplitude (see Fig.l) are regularized by phase space mappings |19[ to ensure 
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an efficient application of VEGAS integrator |20 

While in the Standard Model the cross section of hhZ production |21 
of order 2 • 10 _1 fb at the Higgs boson mass 100 GeV and slowly decreas- 
ing when the mass of Higgs boson increases, the picture in the two-doublet 
MSSM sector is strongly changed by the availability of resonant production 
mechanisms, when the decays of on-shell H — > hh and A — > Zh become 
possible. We show the dependence of total cross sections in the channels 
(18) from the masses of CP-even states rrih and rrin in Fig. 2, 3. In order to 
understand qualitatively the cross section behaviour we show also the h, H 
branching ratio dependence (in the two-body decay channels with the con- 
tribution greater than 1%) from their masses in Fig.4,5. Rapid decrease of 
the total rate at rrih = 60 GeV (m# = 120 GeV) and rapid increase at rrih = 
95 GeV (ma — 190 GeV) are directly connected with the resonant threshold 
of the heavy scalar decay H — > hh (see diagrams in Fig.l). The channel 
e + e~ — > hhZ receives some enhancement at rrih — 95 GeV (m# = 210 GeV) 
when the resonance threshold A — > Zh is opened. 



Our results are qualitatively consistent with the results of JL4], [15], [17 
where somewhat different regions of the two-Higgs-doublet model parameter 
space were explored. Radiatively corrected masses are rather sensitive to the 
input parameter values. At smaller value of tg/3 a mass interval between the 
closing and opening hh thresholds decreases to a few GeV (p^|, tgf3 = 1.5). 

The reactions (18) do not include quartic Higgs boson interaction vertices. 
We calculated the cross-section of the simplest process 

e + e _ — > hhhZ 

(see Fig.2,3), where quartic vertices hhhh and hhhH participate (21 diagrams 
in the unitary gauge). In a very limited region of parameter space the reaction 
has an observable cross-section if the luminosity is high, and the experimental 
reconstruction of multijet events is very efficient. 
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4 Conclusions 



Large increase of the estimate of the possibly achievable integrated luminos- 
ity in the next linear colliders (especially L= 500 fb _1 /year for the TESLA 
project) makes quite realistic the experimental study of Higgs boson self- 
interaction. Such an investigation is especially interesting if the Higgs sector 
of the model includes more than one SU{2) multiplet. Untrivial spectrum of 
scalars leads to the resonant multiple Higgs boson production mechanisms, 
when the final states with 4 or 6 6-jets from their decays will appear with the 
cross sections of one-two orders of magnitude greater than in the SM case of 
only one scalar boson in the Higgs sector. 

For a systematical study of various production channels we derive in a 
compact form a complete set of Feynman rules for the general case of two- 
Higgs doublet model. We demonstrate that in the case of minimal super- 
symmetry, when additional constraints are imposed on the general parame- 
ter space, the interaction vertices are reduced to the well-known vertices of 
the MSSM at the scale M SUS y- Useful connection of LanHEP output in the 
standard lagrangian format of CompHEP input makes possible the following 
efficient calculation of various reactions. 
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Fie 


Ids in the vertex 


Variational derivative of Lagrangian by fields 


H 


H H 


3 M w -Z- S2/3 2( 2s c*-0 2s P+ a vi 2 c a 3 s 2p s M H 2 c s 20 s a 3 M H 2 ) 


H 


H h 


lM w .r w -s 20 *( 12c <*tA-2S a + 4cpi_ilzSp 2M H 2 s 2a s 20 M h 2 s 2a s 2p ) 


H 


H+ H 


m w -L-s^ ( c 3s 2ps a M II ' 2 + c a s 2 ps p 3 M II 2 2s /3+Q /i 2 2 + C{ 3 _ a s 2(i 2 M H ± 2 ) 


H 


A A 


M w -L-s^ ( C /3 3 S2/3S Q A4"i ? 2 + C a S 2( iS(i 3 M H 2 2s p +a ^ 2 + c P - a s 20 2 M A 2 ) 


H 


h h 


2 m w .L-sz * ( M H 2s 2 a s 20 12c Q /i 2 2 s Q + AcptAiSn + 2M h 2 s 2a s 2p ) 


H+ 


11 h 


Myy-JL-s^ ( 2c /3 + aAi?2 C a Cp 3 S 20 M h 2 + S 2 pS a Sp 3 M h 2 + S 2 2 S a - M H ± 2 ) 


A 


A h 


m w -L-s^ ( 2c /3+aA*? 2 c a cp 3 s 20 M h 2 + s 2p s a s p 3 M h 2 + s 20 2 s a _ M A 2 ) 


h 


h h 


3 M w -s w -s 20 ^ ( 2c /3-a 2c /5 + a/ i 12 C a 3 C S 2/3 M h 2 + S 2/ 3S Q 3 S /3 7Vf ft 2 ) 



Table 1. Triple Higgs boson interaction vertices in the general two Higgs 

doublet model 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 


H 


H 


H H 


f Af w 2.^2. S2<3 3 (8s a -p 2 sp +a 2 /2 2 2 - 4(sp ■ c a 3 + cp ■ s a 3 ) 2 s 2 pM H 2 — s 2a 2 s 2 ps a 


-p 2 M h 2 ) 


H 


H 


H h 


(2M H 2 c a 3 s 2 pSp + 2M H 2 cps 2 ps a 3 + M h 2 cp_ a s 2a s 2 p 














H 


H 


H+ IT 


4 M w 2 -s w 2 -s 20 3 ( 4c 2/3 S a-P A«12 + 4S/3+Q H 12 S 2a S 2 p M H 3s 2 pS a -p Sp +a 


M H 2 








— s 2 ps a -p 4 M H 2 — 2c a cp 3 s 2a s 2 ps a _pM h 2 + 2s 2a s 2 ps a s a _pSp 3 M h 2 










-2cp_ a 2 s 2 p 3 M H± 2 ) 




H 


H 


A A° 


4 Af w 2 -s m 2 -s 2/3 3 ( 4c 2/3 Sa -P ^12 + 4S/3+Q //12 S 2a S 2 p M H 3s 2 pS a -p Sp +a 


M H 2 








-S 2 pS a -p 4 M H 2 - 2c a Cp 3 S 2a S 2 pS a -pM h 2 + 2s 2a S 2 pS a S a -pSp 3 M h 2 










-2c P - a 2 s 2 p 3 M A 2 ) 




H 


H 


h h 


4 M lv ?-L 2 -s 20 3 (3s 2a 2 s 2 ps a _p 2 M H 2 + s 2a s 2l3 2 M H 2 3s 2a 2 s 2 ps a _p 2 M h 2 










+4cp +a s 2a s 2 ps a -pM h 2 + 2cp- a s 2a s 2 pSp +a M h 2 — 24c a 2 s a 2 n\ 2 










+8cp 2 s p 2 /nl 2 ) 




H 


H+ 


H h 


1 £^ / A 2 o 2 2 3 7i /f 2 

2 M w 2 -s w 2 -s 20 3 (4C2/3S 2Q //i2 2c /3 _ q S 2/3 S a _pfl 12 C/3 S 2a S 2 pS a S a _ p M H 
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Table 2. Quartic Higgs boson interaction vertices in the general two Higgs 

doublet model 
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Table 3. Triple and quartic Higgs boson interaction vertices at the scale 

Msusy 
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Figure captions 

Fig. 1 Feynman diagrams for the process e + e~ — > hhZ 

Fig. 2 Total cross sections for the reactions e + e~ — > hhZ, e + e~ — > 
e + e _ — > u e u e hh and e + e" — > versus the mass of light CP-even Higgs 

boson at y^s =500 GeV 

Fig. 3 Total cross sections for the reactions e + e~ — > hhZ, e + e~ — > /i/iA, 
e + e" — > u e u e hh and e + e~ — > hhhZ versus the mass of heavy CP-even Higgs 
boson at =500 GeV 

Fig. 4 Two-body branching ratios of heavy CP-even Higgs boson 

Fig. 5 Two-body branching ratios of CP-odd Higgs boson 
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Figure 4: 
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